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The eeyy + fir event observed by the CDF at Fermilab is naturally explained 



'-£3 \ by dynamically supersymmetry breaking models and suggests the presence of 

i ' 

Q_i' the light gravitino which can be a warm dark matter. We consider large scale 

structure of the universe in the worm dark matter model and find that the 



warm dark matter plays almost the same role in the formation of the large 

h : 

scale structure as a cold dark matter if its mass is about 0.5keV. We also 

study the Ly a absorption systems which are presumed to be galaxies at high 

redshifts and show that the baryon density in the damped Ly a absorption 

systems predicted by the warm dark matter model is quite consistent with 

the present observation. 
14.80, 95.35, 98.54.A, 98.65.-r 



1 



Typeset using REVTgX 



Low-energy supersymmetry(SUSY) is a very attractive candidate beyond the standard 



exists the SUSY it must be spontaneously broken. The hidden sector model in N = 1 
supergravity || is widely used for a realization of the SUSY breaking. Although this model 
has many attractive features, it suffers from a serious cosmological problem, i.e. the Polonyi 
problem [|J. There is, however, no such a problem in an alternative model H where the 
SUSY is broken dynamically by some new strong gauge interactions. In this class of models 
the SUSY breaking is mediated to the ordinary sector by the ordinary gauge interactions 
and another problem in the SUSY standard model, i.e. the flavor changing neutral current 
problem, is also solved automatically. 

This dynamical model predicts the SUSY breaking scale F to be low as (100 — 1000)TeV 
and the gravitino mass m 3 / 2 in the range of lOeV — IkeV. Therefore, the usual lightest 
SUSY particle (LSP) decays into the gravitino [] and cannot be a stable cold dark matter 
(CDM) in the universe. Instead of it the gravitino is a true LSP and can form a dark matter. 
Because its mass is so small, the gravitino has larger a velocity dispersion than the CDM. 
Such a type of dark matter is called a warm dark matter (WDM) ||. 

It has been, recently, pointed out [0] that the dynamical SUSY-breaking model naturally 
explains the eeyy + fa event observed by the CDF experiment |§ . The event is explained 
by sequent decays |7j]; e~(e + ) — > e _ (e + ) + B and B — > 7 + G where e, B and G are selectron, 
bino and gravitino, respectively. The decay length of the bino into a photon is given by 



where is the bino mass which should be (38 — 100)GeV 0. For m 3 / 2 < 0.5keV and 
Mg = lOOGeV, the decay length is less than about 5m, which is consistent with the bino 



The longitudinal component of the gravitino (Goldstino) couples to matter with strength pro- 



portional to F 1 ~ 1/ jm^nMp (M p : Planck mass). Thus the decay is very fast for the light 



model, since it provides a natural solution to the gauge hierarchy problem If there 




(1) 




gravitino. 
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decay inside the CDF detector. Thus, the light gravitino is well motivated. 

In this letter, we show that the warm (gravitino) dark matter whose mass is about 
0.5keV 2 plays almost the same role in the large-scale structure formation as CDM. We 
also study damped Ly a absorption systems which are presumed to be the progenitors of 
present-day spiral galaxies in the WDM model and find that the observed mass of neutral 
hydrogens in Ly a systems is consistent with the prediction by the WDM model. Colombi 



et al. |L0| has recently studied the large scale structure formations by WDM. However they 
have considered a very light WDM (~ lOOeV). Such a light WDM may cause a serious 
problem against the damped Ly a systems. 

Before the universe becomes colder than the gravitino mass (T > 10 7 (m3/2/lkeV)K), 
the gravitino behaves as a relativistic particle. Therefore the free streaming of the gravitino 
smears out small-scale density fluctuations and leads to a sharp cutoff in the power spectrum 
of the density fluctuations. The cutoff scale ( = free streaming scale) is given by 

4 = 0.2 ^- (fitf^Mpc, (2) 



uoo, 

where g is the effective number of particle degrees of freedom when the gravitino decou- 
pled (g ~ 200 for the particle content of the minimal SUSY standard model and here- 
after we take g = 200), h the Hubble constant in units of lOOkm/s/Mpc and Q the 
present density parameter of the gravitino which is related to the gravitino mass 7B3/2 
by Qh 2 = (g/100) _1 (m 3/ /2/keV) . Since we only consider a gravitino dominated universe, 
Q ~ Qq where Qq is the total density parameter at present. Assuming a scale invariant 



Harrison-Zeldovich spectrum, we can write the power spectrum P(k) of WDM as [|TT 



P{k) = Ak\T(k)\ 2 , (3) 



2 If we take, e.g. m.3/2 — 0.3keV, we get cr^ ~ 2m. In this case the density parameter of the bino 
is ~ 0.5 and hence we need an additional contribution to O, to get the flat universe. However the 
results in the text are unchanged if the additional contribution comes from CDM || . 
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T{k) = exp 



kR is (kRf. 
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ToW, (4) 



^Ml±^) [1 + 3 , 99 (5) 

+ (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ]- 1/4 , (6) 



where g is defined as q = k/Vl^h 2 / exp(— fi^ — y /7/0.5f2£/fio)Mpc _1 taking into account the 
dependence on the baryon density fig [0, T(k) is the transfer functions for WDM, and A 
the normalization constant which is determined by COBE DMR 4 year data |TB| . Notice that 
for a CDM-dominated universe the power spectrum Pcdm is given by Pcdm(^) = Ak\T (k)\ 2 . 
The WDM power spectrum for = 1, h = 0.5 are shown in Fig. 1 together with the CDM 
power spectrum with the same cosmological parameters. Here, we have taken 7713/2 = 0.5keV 
corresponding to the Qq — 1 universe. Since the cutoff scale is relatively small (~ 0.3Mpc) 
the power spectrum relevant for the large-scale structure (k ^ 1/iMpc^ 1 ) is almost the same 
as the CDM one. (This contrasts the hot dark matter spectrum which has a cutoff of the 
order of O.lMpc -1 .) Therefore, the WDM model with Qq = 1 has the same problem as 
CDM one, i.e. the shape and the magnitude do not fit the observational data from the 
galaxy surveys which are also shown in Fig. 1. In particular, the amplitude of the power 
spectrum normalized by COBE is too large for k = (0.03 — 0.3)/7Mpc _1 . The amplitude 
contradicts not only the galaxy surveys but also the recent analysis of velocity fields |15| 



(shaded region) which could directly reflect the mass distribution. If we take smaller value 
for A, the power spectrum P(k) better fits to the data. Since the tensor mode fl6| 



or 



isocurvature mode [0 may significantly contribute to ST/T in COBE scales, it is possible 
that the actual normalization of A is smaller. For example, we show the power spectrum 
normalized by erg = 0.8 with the same cosmological parameters in Fig. 1. Here as is the 
mass overdensity within spheres of radius 8/i _1 Mpc. Notice that the COBE normalization 
gives <7s = 1.2. As is seen in Fig. 1 P{k) with cr 8 = 0.8 is in a good agreement with the 
velocity field data. 

The difference between WDM and CDM is more significant for galaxies or smaller sys- 



4 



terns. The existence of the cutoff in the WDM spectrum delays the galaxy formation com- 
pared with in the CDM case. Damped Ly a absorption systems observed in QSO spectra are 
important since they give us information about galactic systems in the early universe. It is 
presumed that the damped Ly a absorptions observed in QSO spectra are due to neutral hy- 
drogens contained in galactic systems at high redshifts (z ~ 1 — 4). Therefore, the observed 
damped Ly a absorptions can give us interesting information about baryons contained in 
the galactic systems and can set a constraint on galaxy formation models. In fact, this 
constraint is very stringent for the mixed dark matter (MDM = hot + cold dark matter) 



model |19| , |20| since few galaxies are formed at high redshifts in the MDM model. 

Here we study the damped Ly a constraint on the WDM model. Following ref. ||19|| , we 
use the Press-Schechter theory |2l] to estimate the comoving number density N(z, M) of 
the dark matter halos with mass between M and M + dM at redshift z: 



N(z, M)dM 



'2 po S c 



ttMD^z) 
x exp 



1 da{M) 
a 2 (M) dM 



(7) 



2a 2 {M)Dl{z 

where po is the mean comoving mass density, S c is the overdensity threshold for the collapse 
( = 1.68, corresponding to the prediction of the spherical collapse model), -Di(z) is the 
function for the growth of the perturbations (D\(z) = (1 + z)~ l for Qq — 1) and a 2 (M) is 
the rms mass fluctuation in a top-hat window with radius r M = [M/(47rp /3)] 1//3 , given by 

1 



a\M) 



2tt 2 
3 



P(k)W 2 (kr M )k 2 dk, 



(8) 



W(x) = ^-(sinx — x cosx). 



Since we assume that the damped Ly a absorptions are due to the neutral hydrogens in 
galactic systems, we need identify the halos with a certain mass range as galaxies. For this 
purpose, it is convenient to use circular velocity v c which is related to M by 



M = 2.45 x 10 n M Q (l + z)~ 3/2 (n h 2 )- 1/2 
x(w c /100km/s) 3 ^- 3 . 



(9) 



For the spiral galaxies the circular velocity v c is in the range of 100 — 250km/s. Smaller 
object may also contribute to the absorption. However the lower limit may not be less 
than 50km/s since baryons in such small halo cannot be cooled enough to form a gaseous 
disk [fm . The upper limit is also uncertain since gaseous disks might survive when halo 
merging occurs |19|| . In Fig. 2(a) we show the density parameter of baryons in galactic 



systems (fin) predicted in the WDM model for Q Q = 1, h = 0.5, and Qb = 0.05 together 
with observational data Notice that the adopted Qb is consistent with the prediction 
by big bang nucleosynthesis (tt B h 2 = 0.0125 ± 0.0025 H). In this case, the predicted Q,£> 
with v c = 100 — 250km/s is above the data. Since some fraction of baryons become stars 
and do not contribute to the absorption, the data should be taken to be lower limit to 
Qo- Therefore the prediction by the WDM model is quite consistent with the observation. 
For comparison, we also show the predictions by CDM and MDM in Fig. 2(a). The MDM 
model (i\ot — 0.3 and Q co \d — 0.7) can explain the data of the large scale structure in the 
universe better than the CDM model |23| . However, as already mentioned, since the galaxy 



formation in the MDM model seems too late, it is difficult to explain the damped Ly a 
absorption [T^f2^]. One can also consider the (warm + hot) dark matter model which 



may explain the large scale structure of the universe. However, this model has the same 
difficulty in explaining the damped Ly a absorption systems as the MDM model. In the 
case of CDM, the predicted Qo is larger than that predicted by WDM at higher redshifts. 
Thus, the survey of the damped Ly a systems at z > 5 may possibly distinguish two models. 

As is seen before, the WDM power spectrum with Q = 1 and COBE normalization 
gives poor fit to the data at large scales, which leads us to consider the WDM power 
spectrum with smaller normalization constant. In Fig. 2(b) the predicted Qd is shown for 
Q = 1, h = 0.5, Qb — 0.05, ag = 0.8. In the figure, it is seen that the WDM model is still 
quite consistent with the observation even if we take the smaller normalization for the power 
spectrum. 

If the light fermion forms dark matter in a galactic halo, its phase space density in the 



galactic core might be larger than that allowed by the Fermi statistics [p6 |. This phase space 



constraint puts a stringent constraint on the mass of the dark matter fermion. From the 
study on the stellar motions in dwarf galaxies, the mass should be larger than about 500eV. 
Since the analysis of dwarf galaxies may contain systematic errors, this constraint may not 
be taken seriously. In our case, the dark matter is gravitino and its mass is about 500eV for 
Q = 1 and h = 0.5. Therefore the phase space constraint is satisfied even if we take the 
constraint from the dwarf galaxies. 

In summary, the e&yy + fa event observed by CDF is naturally explained by dynamical 
SUSY-breaking models and suggests the light gravitino. The light gravitino is a LSP and 
can be a WDM. We study the formation of the large scale structure of the universe and the 
adamped Ly a absorption systems in the WDM model. It is found that the WDM plays 
almost the same role in the formation of the large scale structure as the CDM if its mass is 
about 0.5keV. The difference between CDM and WDM becomes more significant when one 
considers the damped Ly a absorption systems which are presumed to be galaxies at high 
redshifts. The baryon density in the damped Ly a absorption systems predicted by the WDM 
model is quite consistent with the observational data. The future observation of damped 
Ly a systems at higher redshifts may distinguish the WDM model from the CDM one. 
Measurements of cosmic microwave background anisotropies may provide another possible 



way to distinguish WDM from CDM |27| although the difference is very small and only 
appears on fine angular scales. We need to wait for new generation satellite experiments(e.g., 
MAP and COBRAS/SAMBA). 
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FIGURES 

FIG. 1. The matter power spectra P(k) in WDM models for the 0,q = 1, h = 0.5 and 0,b = 0.05 
with n = 1 adiabatic fluctuations. The power spectra normalized by the COBE 4 year data(solid) 
and by as = 0.8 (dotted) are plotted. An adiabatic CDM model with COBE 4 year normalization 
is also plotted (dashed). The observational data of galaxy surveys are taken from Peacock and 
Dodds jn[]. Shaded regions are the best fitted value of Mark III catalog of peculiar velocities of 
galaxies by Zaroubi et al. [jll| with 30% errors. 

FIG. 2. (a)Evolution of the Density parameter of the baryons contained in galactic systems 
whose circular velocity is between 100 and 250km/s using Press-Schechter theory and the power 
spectrum normalized by COBE. Solid, dashed and dotted-dashed curves denote the predictions by 
WDM, CDM and MDM, respectively. Symbols represent for the observational data by Wolfe et 
al. fl8|| .(b) Same as (a) with the power spectrum normalized by a% = 0.8. 
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